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Abstract 

The  problem  of  determining  the  acoustic  field  in  an  inviscid,  isentropic 
fluid  generated  fay  a  solid  faody  whose  surface  executes  prescribed  vibrations 
is  formulated  and  solved  as  a  multiple  scales  perturbation  problem,  using  the 
Mach  number  M  based  on  the  maximum  surface  velocity  as  the  perturbation 
parameter.  Following  the  idea  of  multiple  scales,  new  "slow"  spacial  scales 
are  introduced,  which  are  defined  as  the  usual  physical  spacial  scale 
multiplied  by  powers  of  M.  The  governing  nonlinear  differential  equations 
lead  to  a  sequence  of  linear  problems  for  the  perturbation  coefficient 
functions.  However,  it  is  shown  that  the  higher  order  perturbation  functions 
obtained  in  this  manner  will  dominate  the  lower  order  solutions  unless  their 
dependence  on  the  slow  spacial  scales  is  chosen  in  a  certain  manner.  In 
particular,  it  is  shown  that  the  perturbation  functions  must  satisfy  an 
equation  similar  to  Burgers’  equation,  with  a  slow  spacial  scale  playing  the 
r'^le  of  the  time-like  variable.  The  method  is  illustrated  by  a  simple 
one-dimensional  example,  as  well  as  by  three  different  cases  of  a  vibrating 
sphere.  'The  results  are  compared  with  solutions  obtained  by  purely  numerical 
methods  and  some  insights  provided  by  the  perturbation  approach  are  discussed. 

'The  first  author  was  .supported  by  NASA  Research  Grant  NAG-l-lO!)').  Research  was  also  suppc>rl«d  by  the 
National  Aeronautics  and  Space  Administration  under  NASA  Contract  No.s.  NASi-lHfiO'i  aiul  NASI  lUISO  wliilc 
the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and  F.ngineeriiiR  (U'ASF,),  N  ASA 
Langley  Research  Center,  Hampton,  VA  1-0001 . 


1.  Introduction 


We  wish  to  describe  the  acoustic  field  in  an  inviscid,  isentropic  fluid 
generated  by  a  solid  body  whose  surface  executes  prescribed  vibrations  When 
the  Mach  number  M,  based  on  the  maximum  surface  velocity  is  small,  it  is 
natural  to  consider  a  perturbation  approach  to  the  problem.  Crow  (1970) 
carefully  analyzed  a  general  class  of  problems  of  this  type,  including 
Lighthill’s  (1952,  1954)  acoustical  analogy  approach,  involving  quadrupole 
sources,  and  Ribner’s  (1962)  approach,  based  on  monopole  sources.  He  showed 
how  this  class  of  problems  could  be  treated  as  a  singular  perturbation  problem 
by  the  method  of  inner  and  outer  asymptotic  expansions.  Whitham  (1956,  1974) 
examined  this  problem  in  the  high  frequency  limit  by  using  a  perturbation 
approach  similar  to  Lighthill’s  (1949,  1961)  method  of  strained  coordinates. 

He  essentially  "strained"  the  time  coordinate  and  obtained  expressions  for  the 
acoustic  field  valid  in  the  far  field  for  high  frequency  excitations.  The 
approach  we  shall  use  is  not  restricted  by  the  high  frequency  assumption, 
although  our  results  do  reduce  to  Whitham’s  results  (to  leading  order  in  M)  in 
the  far  field  and  at  high  frequencies.  In  addition,  the  analysis  we  shall 
present,  being  based  on  multiple  spacial  scales,  provides  insight  into  the 
effects  of  the  nonlinearity  on  the  form  of  the  solution,  as  well  as  insights 
irto  the  interaction  of  the  value  of  the  Mach  number  and  the  angular 
dependence  of  the  solution  (see  section  9). 

In  section  2  we  formulate  the  mathematical  problem  we  wish  to  solve  and 
then  outline  the  multiple  scales  approach  in  section  3.  We  apply  our  method 
to  a  simple  one-dimensional  example  in  section  4.  In  section  5  we  formulate 
the  problem  for  a  vibrating  sphere  and  then  determine  the  leading  terms  in  the 
multiple  scales  expansion  in  sections  6  and  7.  In  section  8  we  apply  our 
results  to  three  different  vibrating  spheres  and  then  discuss  our  results  in 
section  9. 

2.  Problem  Formulation 

We  wish  to  describe  the  acoustic  field  in  an  inviscid,  isentropic  fluid 
generated  by  a  solid  body  whose  surface  executes  prescribed  vibrations.  To 
begin,  we  let  the  origin  of  a  Cartesian  coordinate  system  (Xj,X2,X2)  be  fixed 

at  some  convenient  location  inside  the  body.  Then,  in  terms  of  these 


1 


coordinates,  the  equations  of  conservation  of  mass  and  momersturo  are 


(2.1)  dpidt  +  ^•(pu)  =  0, 

(2.2)  dpldx.^  +  p\^.^ldt  +  (u*^;u.j  =  0,  i  =  l,2,3 


In  (2.1)  and  (2.2)  p,  p,  and  t  represent  the  fluid  density,  the  fluid 
velocity  component  in  the  positive  x.  direction,  the  fluid  pressure,  and  the 

time  coordinate,  respectively.  Also,  u  =  u  i  +  u  i,^  +  u  i  and  ^  = 


i.,(d/(9x,  )  +  \Adldx^)  +  i„(5/dx-),  where  i.  is  a  unit  vector  in  the  positive 

Jl  \  2,  £,  o  ^  J 

Xj  direction.  We  also  let  p  =  kp  ,  where  k  and  -7  are  known  constants,  and 
define 


(2.3)  c^(p)  =  dp/dp  =  k7(p)^ 

We  let  pQ,  U,  and  L  be  typical  (constant)  values  for  the  density, 

velocity,  and  length  scales,  respectively,  associated  with  the  flow,  and 
define  nondimens ional  variables  (denoted  by  a  above  the  quantity)  by 


(2.4) 


p  =  p/pQ,  u  =  u/U,  p  =  (p-Pq)/(PqU^) , 
X.  =  x./L  (j=l,2,3),  t  =  t/(L/CQ), 

J  J 


where  =  kp^  and  c^  =  zip^) .  We  now  rewrite  equations  (2.1)  and  (2.2)  in 

terms  of  these  variables  and  then  omit  the  above  the  various  quantities  to 

obtain  the  relations 


(2.5) 


dpidt  +  y^^ipxx)  =  0, 
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(2.6)  dpJdx^  +  +  .M^(u*^)u.j  =  0,  i=l,2.3  , 

where  M  s  is  the  Mach  number  of  the  flow. 

We  now  seek  approximate  solutions  to  equations  (2.5)  and  (2.6),  subject 
to  appropriate  boundary  and  'nitiai  conditions,  which  will  be  formally  valid 
for  small  values  of  M. 


3.  A  multiple  scales  perturbation  solution  of  the  basic  equations 

Following  the  method  of  multiple  scales  (see  e.g.  .Nayfeh)  1973 ) ,  Chapter 

k) 

5),  we  introduce  the  spacial  scales  (variables)  x  ,  k  =  0 ,  1 ,  2 ,  . . .  , 
related  to  x  by 

f  k  1  k“^ 

(3.1)  ’  =  M  X,  k  =  0,  1,  2 . 


-*^(1)  -^(2)  2-^ 

Thus  X  =  X,  X  =  Mx,  X  =  M  X,  etc.  We  then  assume  that  u  and  p  are 


functions  of  these  new  scales,  as  well  as  functions  of  t  and  M,  i.e.  p  = 
nrf  t^O)  v^2)  n  -  -itrf  v^2) 


(k) 

The  method  of  multiple  scales  now  treats  all  of  the  variables  {x^  )  as 

independent  variables  and  seeks  to  determine  p  and  u  as  functions  of  these 

variables.  In  particular,  for  "small"  values  of  M,  we  look  for  solutions  in 
the  form 


QD 


(3.2) 


p  =  1  +  Mp^^^  +  +  ...  =  ,  with  p^^^=  1, 


j=0 


00 


^  -^(0)  -(1)  „2-(2)  ^^(j)uj 

u  =  u  +  Mu  +  M  u  +  ...  =  >  u  •'  M-'  . 

j=0 
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Here  each  of  the  coefficient  functions  ^ '  and  is  independent  of  M,  but 

in  general,  will  depend  upon  t  and  the  spacial  scales  i  e  = 


;^(0)  -(1)  >  ^(j),  -(0)  -(1) 

p  (t,x  ,x  ,...)  and  u''=u-'(t.x  ,x 


}  . 


To  determine  these  coefficient  functions,  we  substitute  the  expansions 
(3.2)  into  equations  (2.5)-(2.6)  and  use  the  relation 


(3.3) 

1  e  ...  . 

where 

We  then  collect  coefficients  of  like  powers  of  M  on  the  left  side  of  each 

equation  and,  hence,  express  the  left  side  of  each  equation  as  a  power  series 

in  M.  We  then  equate  to  zero  the  coefficient  of  each  power  of  M,  since  the 
right  side  of  each  equation  is  zero.  In  this  way,  v.e  obtain  the  following 

system  of  equations  satisfied  by  the  coefficient  functions  ^ ^  and  u^^^: 

(3.4a) 

dp^'^hdt  +  =  0, 

(3.4b) 

du!^^/dt  =0,  i  =  1,  2,  3; 

(3.5a) 

(3.5b) 

dp^^*^hdx[^\  dii^^hdt  =  i  =  1,  2,  3,  k^l . 

(Equations  (3.4)  follow  from  the  terms  in  equations  (2.5)-(2.6)  which  are 

0(M),  while  equations  (3.5)  follow  from  the  terms  in  these  equations  which 

k+l  fk)  f  kl  f  i  i 

are  0(M  ).)  Here  the  functions  F'  and  depend  only  upon  p  •'  with 

i ) 

j<k+l  and  u  •'  with  j<k.  In  particular,  we  find 

(3.6a) 
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(3,6b) 


Op"', -Ox'.'') 


To  solve  equations  (3.4),  we  set 

(3.7)  -^^/dt  and  u-®^=  ■ 

Then  equations  (3.4b)  are  satisfied  for  any  choice  of  ‘{i,  while  equation 
(3.4a)  yields  the  requirement  that  'p  must  satisfy  the  I’lear  wave  equation 

(3.8)  d^p/dt^  -  =  0, 

where  A^^^s  (d/3x|^^)^  +  is  the  usual  Laplacian 

operator  in  the  variables 

From  the  structure  of  equations  (3.5),  we  see  that  these  equations  can 
(in  principle)  be  solved  recursively,  starting  with  k=l.  In  particular, 

equations  (3.5)  are  a  system  of  linear  equations  for  the  unknowns  and 

“^(k) 

u  .  Consequently,  we  can  express  the  solution  to  these  equations  as  the 
superposition  of  a  particular  solution  and  a  homogeneous  solution  to  these 
equations.  The  homogeneous  solution  has  the  same  form  as  the  solution  to 
equations  (3.4)  (see  equations  (3.7)  and  (3.8)).  The  particular  solution 
will,  of  course,  depend  on  the  index  k,  as  well  as  on  the  solutions  for  the 
lower  order  perturbation  coefficients.  In  general,  these  particular  solutions 
have  the  property  (as  we  shall  demonstrate  explicitly  in  the  following 

sections)  that  they  tend  to  decay  more  slowly  in  magnitude  as  |xj->oo  than 

-k(0) 

and  u  Consequently,  if  this  decay  were  to  be  left  unchecked,  the 

perturbation  expansions  (3.2)  would  become  invalid  as  |xl->®.  However,  as  we 

{ lc+1 )  )c ) 

shall  show,  it  is  possible  to  satisfy  the  requirement  that  p  and  u 

should  decay  at  least  as  fast  as  p^^^  and  u^^^  as  |x|-^  by  properly  choosing 

(k+1)  "^(k)  “^(1) 

the  dependence  of  p  and  u  on  the  "slower’'  spacial  scales  x  ,  x  , 

etc.  Thus,  this  method  will  yield  a  perturbation  expansion  (3.2)  which  will 
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be  uniformly  valid  as  jx{-»«!and,  consequently,  will  correctly  represent  the 
far  field  behavior  of  the  acoustical  radiation. 

In  the  following  section  we  shall  demonstrate  some  of  these  ideas  with  a 
simple,  one-dimensional  example,  and  then  proceed  to  a  class  of  three- 
dimensional  problems  in  sections  5-8. 

4 .  A  Simple  One-Dimensional  Example 

As  a  simple  example  to  illustrate  the  general  ideas  of  our  approach,  we 
consider  first  a  flow  which  varies  in  only  one  spacial  coordinate,  such  as 
one-dimensional  flow  in  a  semi-infinite  tube.  We  let  this  one  spacial 
coordinate  be  denoted  by  x^^  =  x,  so  that  equations  (2.5)-(2.6)  become 

(4.1)  dp/dt  +  MO/dx)(pu)  =  0  , 

(4.2)  dpidx  +  p^  +  M^u(^/dx)j  =  0  . 


We  shall  assume  that  u  is  a  specified  function,  say  f(t),  at  x=0  for  all  t^O, 
and  that  ambient  conditions  prevail  at  t=0  for  all  0<x<oc,  i.e. 


(4.3) 


u 


=  f(t)  for  t^O,  and 
x=0 


u  =  0,  p  =  1  ,  for  x>0. 
'  t=0  ' t=0 


Following  the  method  of  multiple  scales  outlined  above,  we  find  it 
convenient  to  define  the  differential  operators  D.  by 

(4.4)  D.  H  d/dx^^\  i  =  0,  1,  ...  . 

We  then  look  for  solutions  for  p  and  u  in  the  form  of  (3.2),  where,  in 

particular,  p^^^  and  u^^^  satisfy  equations  (3.4).  Then  the  general  solution 
to  equation  (3.8)  is  given  by  p  =  F(t-x)  +  G(t+x),  where  F  and  G  are 
arbitrary  functions  of  their  arguments  and  may  also  depend  on  the  spacial 

scales  x^^\  •  Since  we  shall  restrict  our  attention  to  solutions 
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which  propagate  only  in  the  positive  x  direction,  we  set  G  s  0  and  write 


(4.5) 


f ( t-x. 


X 


(2) 


.  ) , 


where  the  exact  form  of  the  function  f  will  be  determined  from  the  boundary 
and  initial  conditions  (4.3)  of  the  problem. 

Using  the  solutions  (4.5),  equations  (3.5)  with  k=l  for  and  u^  ^  ^ 

become 


(4.6) 


ap^^^at  +  DqU^^^  =  -[o^f  +  OQ(f^)] 
au^^^/at  +  DqP^^^  =  -(^D^f  +  ((7-l)/2)DQ(f^)]  . 


The  general  solution  to  (4.6)  for  p''^^  and  u^^^  can  be  written  as 

(4.7) 


In  (4.7),  fj^  is  an  arbitrary  function  of  its  arguments,  which  will  eventually 
be  determined  by  the  initial  and  boundary  conditions  of  the  problem. 

In  order  for  p^^^  and  to  grow  no  faster  that  p^^^  and  u^^\ 

respectively,  as  x  becomes  large,  we  must  require  that  the  term  in  brackets  in 
equations  (4.7)  vanishes,  i.e. 


p^^^  =  -x[D^f  +  ((l+7)/4)DQ(f^)]  +  fj(t-x,  x^^\  x^^\  ...), 

=  -x[D^f  +  ((l+'y)/4)DQ(f^)) 

+  ((7-3)/4)(f^)  +  fj(t-x,  x^^\  ...). 


(4.8) 


Djf  +  ((l+7)/4)DQ(f^)  =  0  . 


Condition  (4.8) 


is  an  equation  to  determine  the 


behavior 


of  f  and. 
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hence,  the  behavior  of  both  and  Equation  (4.8j  can  be 

expressed  in  terms  of  Burgers’  equation,  with  ( { 1 +'yi /4  ,)x' ^  '  playing  the  role 
of  the  "time"  variable  and  x  playing  the  role  of  the  spacial  variable.  In 
particular,  the  solution  to  (4.8)  which  satisfies  the  "initial"  condition  that 

f  =  f(t-x)  when  x^^^=0  can  be  expressed  as 


(4.9)  f  =  f(T),  where  r  =  t-x  +  ( (7+1 ) /2)x* ^ ^f (r  )  . 

As  an  application  of  these  results,  we  let  the  (dimens i ona 1 j  fluid 
velocity  at  x=0  be  u  =  Usin(iJt).  We  then  use  U  as  the  typical  velocity  of  the 
flow  and  let  L  =  Cq/w  be  the  typical  length  associated  with  the  flow.  (Note 

that  k  =  1/L  =  w/Cq  is  the  linear  wave  number  for  this  example.;  Then  the 

boundary  condition  for  the  flow,  when  expressed  in  terras  of  the  nond iraens i ona  1 
variables  (2.4),  becomes 


(4.10)  u  =  sin(t),  at  x=0  ,  for  all  t^O  (boundary  condition). 


Using  the  expansions  (3.2)  in  (4.10)  we  find  that  the  individual  coefficients 
and  satisfy  the  initial  and  boundary  conditions 


(4.11) 


(4.12) 


0 

0 


(0)  .  , 

u  =  Sint 


0 


at  t=0  (for  all  x^O)  for  k  =  0,  1,  2, 


at  x=0  (for  all  t>0)  for  k  =  1,  2,  3, 


From  equations  (4.5),  we  find 


(4.13) 


f(t-x. 


.  .  .). 
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If  we  were  to  terminate  our  perturbation  expansion  at  this  point  and  set  x 

for  i>l  equal  to  zero,  then  the  conditions  ( 4 . 1 1 j-( 4 . 12)  would  yield  the 
results 


(4.14) 


(1) 


sin(t-x),  for  0^x;^t  , 

0  ,  for  x>t  . 


This,  of  course,  is  just  the  classical  wave  solution  of  linear  acoustics. 
However,  if  we  continue  the  perturbation  solution  outlined  above  and  use  the 
solution  (4.9),  we  find 


(4.16) 


-  ((3-7)/4)sin^(t-x), 

=  ( (7-3) /4)  l^sin^(r)  -  sin^(t-x)j, 


for  0^x;^t,  with  p^^^  =  s  0  for  t>x. 

Combining  the  expressions  we  have  obtained  for  p^^^  and  p^'^\  we  can 

write 

(4.17)  p  =  1  +  Msin(r)  +  M^( (3-7) /4)sin^( t-x)  +  O(M^)  +  O(M^x), 

u  =  sin(r)  +  M( (7-3) /4) ^sin^(r )  -  sin^(t-x)j 
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+  +  OC.M'^xJ, 

for  Oi^x^t,  and  p=l  and  u=0  for  x>t.  Here  r  is  defined  by  equation  (4.1,5bj. 

To  illustrate  these  results,  in  Figure  1  we  have  plotted  approximations 
to  (p-1)  as  a  function  of  distance  from  the  source  at  a  fixed  time.  The 
multiple  scales  approximations  (solid  lines)  were  obtained  using  the  first 
three  terms  on  the  right  side  of  the  first  equation  in  (4.17)  and  are  plotted 
for  two  different  values  of  M.  In  Fig.  1(a),  M=0.01  and  the  resulting  wave 
lies  close  to  linear  solution  (4.14)  (short  dashed  lines),  with  only  a  small 
distortion  due  to  the  nonlinearities  in  the  governing  equations  being  visible. 
In  Fig.  Kb),  M  has  been  increased  to  0.04.  For  this  case,  the  gradual 
steepening  of  the  wave  as  x  increases  is  evident,  with  the  wave  apparently 
approaching  an  "N-wave"  as  x  increases.  An  approximation  based  on  Wliitham’s 
(1974)  first  order  approximation,  obtained  using  essentially  the  method  of 
strained  coordinates,  is  very  close  to  (but  not  identical  to)  the  multiple 
scales  solution,  and  hence  is  not  shown  on  these  plots.  In  these  figu*-  ;  we 
have  also  plotted  a  solution  (circles)  to  the  basic  equations  (4.1)-(4.2) 
obtained  by  purely  numerical  means.  This  numerical  solution  was  obtained  by 
writing  the  basic  equations  in  characteristic  form  and  then  using  a  MacCormack 
predictor-corrector  technique  on  the  resulting  equations.  (We  thank  Dr. 

Willie  R.  Watson  of  the  NASA  Langley  Research  Center  for  carrying  out  these 
calculations  for  us.)  As  the  figures  illustrate,  there  is  very  good  agreement 
between  the  multiple  scales  solution  and  th^  numerical  solution  for  these 
values  of  the  Mach  number.  We  shall  comment  further  on  these  results  in 
section  9. 

5.  Three-Dimensional  Flow  -  The  Vibrating  Sphere 

We  now  wish  to  describe  the  acoustic  field  in  an  inviscid,  isentropic 
fluid  outside  a  spherical  region,  when  the  normal  component  of  the  fluid 
velocity  is  specified  on  the  surface  of  this  region.  It  is  convenient  to 
think  of  this  spherical  region  in  one  of  two  ways.  First,  we  may  regard  it  as 
an  actual  sphere,  whose  surface  executes  prescribed  vibrations,  with  the 
boundary  conditions  applied  on  the  time  averaged  surface  of  the  sphere.  Under 
this  interpretation,  the  results  which  follow  are  extensions  of  several 
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"classical"  problems  in  acoustics.  Ai terna t i ve i y ,  we  may  regard  the  region  as 

a  ^mathematical )  sphere  which  encloses  a  bounded  region  of  complicated  fluid 

flow.  With  this  interpretation,  our  results  may  be  applicable  to  the  problem 

of  determining  the  acoustic  field  radiated  from  bounded  (perhaps  lurhuienti 

fluid  flows.  In  particular,  our  results  will  show  that  only  the  normal 

component  of  the  flow  velocity  needs  to  be  specif itu  on  this  spherical  surface 

in  order  to  determine  the  radiated  acoustical  field 

To  begin,  we  let  the  origin  of  a  ( nond  imens  i  ona  1  ;  Cartt-sian  coordinate 

system  (x  ,x  ,x,)  coincide  with  the  average  position  of  the  center  of  the 
X  ^  o 

sphere,  whose  average  radius  is  a  constant  a.  We  then  introduce  .spherical 
coordinates  (see  Figure  2(a))  related  to  by  the  relations 

Xj=  rs  in{^}cos{^) ,  rs  in(^)s  in(^) ,  x^=  rco.sS^l.  IHtre  all  lengths  have 

been  nondimens i ona 1 i zed  by  referring  them  to  the  average  radius  a  of  the 
spherical  region.)  Then,  in  terms  of  these  coordinates,  equations  ^2.5)-(2.6) 
become 

(5.1)  dp/dt  -*■  m||(  1  / r^)  (5/dr )  (  r'^pu^ )  +  ( 1  / rs inf  0  )  i (  5/5^ )  (  s  i n<  0 .)pu^ i 

(  1  /  r  s  i  n  (  ^ ) )  f  5/5r  )  ( pu  , )  j  =  0 , 

V 

(5.2)  dp! dr  +  p^  +  M^ju^(5u^/5r )  +  ( u^/ r) (5u^ /55 ) 

+  (u^/rsin(0)  )(5u^/5‘|)  -  (u^)^/r  -  (u^)^/rj|  = 

(5.3)  dptdB  +  rp^~^|H5u^/5t  +  M^|^u^{5u^/5r)  +  (u^/r )  (du^/d^  i 

O  '  [ 

+  (u^/rsin(d) )(5u^/5^)  +  (u^u^i/r  -  ( u^)  co t ( 0 ) /r j  |  =  0, 

/ 

(5.4)  dpld^f  +  (rsin(fl))p^  ^|M5u^/5t  +  M^j^u^(5u,^/5r)  +  (u^/r)(5u^/5fl) 

+  (u^/rs in{ ^) ) (5u^/5^)  +  {u^u^)/r  +  (u^u^)cot(0) /r 


In  ( 5 . 1 )  *-{  5. 4 )  ,  u  ,  Urt,  u,  represent  the  fluid  Vflocitv  componeril  s  in  the 
T  tf  Y 

positive  r,  0  and  directions,  respectively 

We  shall  assume  that  is  a  specified  function,  say  t).  at  r=l  fur 

all  t^O,  and  that  ,  u^,  u^  and  p-1  are  all  specified  to  be  zero  at  t=0  for 

al 1  l<r<®,  i . e . 


(5.5) 


u  =  Vid.^.tJ.  (u  u^.u, ) 

r=l  r  (7  y 


=  (0,0,0)  and  p  =  1 
t=0  t=0 

r>l  r>l 


Following  the  method  of  multiple  scales,  we  introduce  the  spatial  scales 
fp  i  =  0,  1,  2,  ...  ,  related  to  r  by 


(5.6) 


ry  =  r.  r.  =  M  r.  1  =  1.  2, 


Then,  for  "small"  values  of  M,  we  look  for  solutions  for  the  density  and 
velocity  components  in  the  form 


(5.7) 


p  =  1  +  +  .  .  .  =  ^  ,  with  1 

j=0 


(0)  „  (1)  u2  (2) 

u=u  +Mu  +Mu  +... 
r  r  r  r 


j=0 


with  analogous  expressions  holding  for  u^  and  u^.  Here  each  of  the 
coefficient  functions  p^'^^  ’  independent  of  M,  but,  in  general 


will  depend  upon  t  and  the  spacial  scales  rp  i.e.  p  ^  = 

p^-^^(t,ry,rj^,r2,...),  uj^-*  u|^j  ^  ( t ,  r^ ,  rp  Cg,  .  .  . ) ,  ...  . 

To  determine  these  coefficient  functions,  we  substitute  the  expansions 
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(5.7)  into  equations  (5.1)-(5.4)  and  use  the  relation 


(5.8)  dldr  =  D„  +  MD.  +  M^Dr,  +  ...  ,  where  D.  ^  di  dr . 

0  12’  j  J 

We  then  collect  coefficients  of  like  powers  of  M  on  the  left  side  of  each 
equation,  as  described  in  section  3,  and  find  the  following  system  of 

equations  satisfied  by  ^  ^ ' 


(Equations  (5.9)  follow  from  the  terms  in  equations  (5.1)-(5.4)  which  are 
0(M),  while  equations  (5.10)  follow  from  the  terms  in  these  equations  which 

are  0(M^^^).)  Here  the  functions  and  G^^^  depend  only  upon 

with  j<k+l  and  u^^\ 

In  the  next  sections  we  shall  show  how  equations  (5.9)-(5.10)  can  be 
solved  recursively. 
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6.  Solution  for  the  lowest  order  perturbation  coefficients 

The  form  of  equations  (5.9)  suggests  that  we  look  for  a  solution  for  the 
lowest  order  perturbation  coefficient  functions  in  the  form 

(6.1)  =  -dipidi,  d-pldr, 

{\lT)dpld9 ,  (l/Ts\n{6))dtp/d4', 

where  p  is  a  smooth  function  of  its  arguments,  which  still  needs  to  be 
determined.  Using  equations  (6.1),  we  see  that  equations  (5.9b)-(5.9d)  are 
satisfied  for  any  choice  of  p,  while  equation  (5.9a)  leads  to  the  requirement 
that  tp  must  satisfy 

(6.2)  ^(pldt^  =  {llr‘^){dldr)^^dipldr^ 

+  (.llT^sin{6)){dldB){sin{$)d<pldd^  +  (l/r^sin^(^))d^p/^^, 

which  is  just  the  usual  (linear)  wave  equation.  (Here  we  have  used  the 
relation  ~  dldr.)  This  equation  has  solutions  of  the  form 

(6.3)  p  =  F^(r,  t)F^(cos(5)  )cos(m^-/?) ,  for  n  =  0,  1,  2,  ...  , 

and  m=0,  1,  ...  ,  n, 

where  is  the  associated  Legendre  polynomial,  /?  is  an  arbitrary  constant, 
and  satisfies  the  differential  equation 

a^F  Idi^  =  (l/r^)(a/ar)fr^aF  /dr]  -  (l/r^)n(n+l)F  . 
n  n  ;  n 

The  outgoing  solution  for  F^  can  be  expressed  as 
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(6.4) 


n+1 


F  (r,t)  =7^8  •'^(t-r+D/r-', 

n  / . -  -  n ,  j 


J=1 


%,!  =  -I’  %,j+l  '  (<n*l-i)(i«)/(2j))a^_.,  i  =  1,  2, 

where  f  is  an  arbitrary  (smooth)  function  of  its  argument. 

If  we  now  apply  the  boundary  condition  (5.5) 


.  n, 


00 


we  see  that  we  can  set 


r=l 


=  D 


r^,  r^, 


j=0 


r=l 


u 


.(0) 


(6.5) 


'r=l 


=  t) ,  for  j=0, 


.(j) 


u.'  =  0  ,  for  j^l  . 

r=l 


We  shall  assume  that  V  can  be  expressed  as  (a  linear  combination  of  terms  of 
the  form) 


(6.6) 


=  Q(  t)P^(cos(^)  )cos(m^-/?) , 


where  Q(t)  is  a  specified  function  of  t,  and  n  and  m  are  non-negative 
integers,  with  O^m^n.  Using  this  expression  for  V,  along  with  the  condition 

(6.5)  on  as  well  as  the  expressions  (6.3)  for  <p  and  (6.1)  for  we 

find  that  f  must  satisfy  the  condition 

n 

(6.7)  f^'^'^^^t)  -  ^[(j(j+l)+n(n+l))/(2j)]a^  ^  ( t) 

j=l 
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~  (n+1 )a  . f ( t )  =  Q( t )  . 

n,n+l 

Equation  (6.7)  is  an  inhomogeneous,  linear  ordinary  differential  equation  of 
order  n+1,  with  constant  coefficients,  for  the  unknown  function  f(t).  For  a 
few  small  values  of  n  this  equation  becomes 

n=0:  f’+  f  =  Q(t); 

n=l;  f”+  2f’  +  2f  =  Q(t); 

(6.8)  n=2:  f  ’  ”+  4f  ’  +  9f ’+  9f  =  Q( t) ; 

n=3;  f””+  71’”+  27f”+60f’+  60f  =  Q(t); 

n=4:  f’””+  llf””+  65f”’+  240f”+  525f’+  525f  =  Q(t). 

Once  f  has  been  determined,  we  can  express  the  solution  for  ^  as 

'  n+1 

(6.9)  y  ' a  .f  ^(t’ )/r^  p”(cos(d)  )cos(m\l'-^) , 

n, j  n 

^  j=l 

where  t’  =  t-r+l.  In  general,  f  will  also  depend  on  the  "slow"  spacial  scales 
r^,  Tg,  although  this  dependence  has  not  been  denoted  explicitly. 

As  an  illustration  of  these  results,  we  consider  the  special  case  when 
Q(t)  =  sin(wt).  Then  the  steady  state  (periodic)  solution  to  (6.7)  (which 
can  be  viewed  mathematically  as  a  particular  solution  to  this  equation)  can  be 
expressed  as 

(6.10)  f(t)  =  Asin(wt-a),  when  Q(t)  =  sin(wt), 
where,  for  a  few  small  values  of  n,  A  and  a  are  given  by 

n=0:  A  =  {X-HtP")  a  =  tan  ^(w), 

n=l;  A  =  (4+w^)  a  =  tan”^  ^2u;/(2-w^) j  , 
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n=2:  A  =  (81+9w^-2w‘*+w®)  a  =  tan"”^  |^w(9-ij^) / ( 9-4w^)j  , 

n=3:  A  =  (3600+360u;^+9w'^-5w®+w®)“^^^, 

a  =  tan"^  j^a;(60-7w^)/(u'*-27u;^+60)j  , 

n=4:  A  =  (275625+23625w^+900w^-5w®-9«®+u;^®)''^^^, 

a  =  tan"^|^w(525-65a;^+u;'^)/(525-240a;^+lla;^)j  . 


Solutions  to  the  homogeneous  version  of  equation  (6.7)  are  of  the  form  f(t)  = 

e^^,  where  the  real  part  of  A  is  negative.  Hence  the  corresponding  solutions 
for  f(t)  are  not  periodic  in  t  and,  in  fact,  they  all  decay  to  zero  as  time 
increases . 

Before  we  consider  the  higher  order  perturbation  coefficients,  we  can  use 
the  results  above  to  calculate  a  first  approximation  to  the  radial  component  I 
of  the  sound  intensity,  defined  by 

(6.11)  I  =  PqCqU<Mp*u^>, 

where  the  symbol  "<g>"  denotes  the  time  averaged  value  of  the  quantity  g. 

Using  the  expressions  (6.1),  (6.3)  and  (6.4)  above  we  find,  for  the  case  when 
Q  is  given  as  in  (6.10),  that 

2 

I  =  pQCQU^w^’°^^A^|^P^(cos(fl))cos(mf-/3)j  /2r^  +  0(l/r^),  as  r-»«. 


For  w  =  ka  «  1,  we  find  that  the  maximum  intensity  of  I,  which  we  denote  by 
I^,  for  different  values  of  n,  is  given  by 


n=0;  Iq  =  [pQCQU^(ka)^/2](a/r)^  +  Oiika)"^/?^)  +  0(l/r^); 

n=l:  =  (^pQCQU^(ka)^/2](a/r)^((ka)^/2^]  +  0((ka)®/?^)  +  0(l/r^); 
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n=2:  =  [/?QCQU^(ka)^/2]{a/r)^[(ka)'^/3^]  +0({ka)®/?^)  +0(l/r^). 


In  general,  for  any  positive  integer  n,  the  ratio  of  the  maximum  intensity  I 

n 

to  the  maximum  intensity  of  the  pulsating  sphere  approaches  (ka)^"/{n+l)^, 
as  ka  approaches  zero. 


7.  Solutions  for  the  higher  order  perturbation  coefficients 
Equations  (5.10)  with  k=l  can  be  written  as 


(7.1a) 


(7.1b) 


Idi  +  (l/r^)DQ|^rV^^] 

+  (l/rsin(d))0/atf)(sin(d)u^^^)  +  (l/rsin(5)  = 

-[(l/r^)DQ(r^/3^^^u[®^)  +  Dju[®^ 

+  (l/rsin(d))(d/d(?)(sin(^)p^^^u^®^) 


+  (l/rsin(d))(d/d^)(p^^^u^®^], 


+  du^^^dt  = 


+  (u^^^/r)(du^^^/5fl)  +  (u^*^^ /rsin(0) ) (^1^^^ /d^) 


(7.1c) 


(7. Id) 


-  (u^®^2/r  -  (u^®^)2/rl. 


dp^^^ Id6  +  rdu^^^/5t  = 


-r[i2-,)p^^W;hdt  . 


+  (u^®^r)(du^®^/a0)  +  (u^°^rsin((l))(du^°Vai^) 


+  (u[°^u^®^/r  -  (u^°^)^cot(fl)/r]  , 


tdi’  +  rsin(5)^^^^/dt  = 
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-(rsin(0)) 

+  (uj.^^u^‘^^)/r  +  (u^®^u^®^cot(^)/r). 


To  begin  our  construction  of  solutions  to  equations  (7.1),  we  first 
examine  the  behavior  of  the  right  sides  of  these  equations  as  r-^.  Using 
equations  (6.1),  (6.3),  and  (6.4)  we  find 

=  -dtpidi  =  (l/r)f^"'''^^P^(cos(fl))cos(m^-/3)  +  0(l/r^), 

(7.2)  =  dpIdT  =  (l/r)f^"^^V(cos(<9))cos(mM)  +0(l/r^), 

=  {llr)dplde  =  (l/r^)f^"^P“’(cos(^))(sin(^))cos(m\;'-/?)  +  0(l/r^), 

u^^^  =  (l/rsin(^) )3p/3^ 

=  (m/r^sin(^))f^"^P^(cos(fi))sin(in\;'-/3)  +0(l/r^),  as  r-H». 

(Here  it  is  understood  that  the  argument  of  f  is  t’=t-r+l.)  Then,  using  these 
expressions,  we  can  write  equations  (7.1)  as 

(7.3a)  dp^^hdt  +  (l/r^)DQ(r\J^^] 

+  (l/rsin(d))(5/3^)(sin(0)u^^^)  +  (l/rsin(^)  )3u^^^/d\f'  = 

-D^uJ.®^  +  (2/r^)f^"''^^f^""'^^[pJJ(cos(0))cos(rai^-^))^  +0(l/r^), 
(7.3b)  = 

2 

+  ((7-l)/r^)f^'^'^^^f^"''^^(p^(cos(ff))cos(m^^-^))  +0(l/r^), 

(7.3c)  dp^^hde  +  = 
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{(l-7)/r^)  P™(cos(5))P^’(cos(^))sin(^)cos^(m^-/?)  +  0(l/r^), 

(7.3d)  +  rsin(5)du^^^dt  = 

2 

m(  (7-I ) /r^)  ^P®(  cos(5 ) )  j  cos(in^-/?)sin(mv!'-/3)  +0(l/r^). 

We  now  examine  the  effect  on  the  solutions  each  of 

the  terms  on  the  right  sides  of  equations  (7.3).  In  particular,  we  must  first 
examine  the  behavior  of  the  terms  and  D,u*^^^  as  r->CD.  Now.  if  we  assume 

that  and  are  0(l/r)  as  r-^,  then  these  terms  will  give  rise  to 

terms  in  p^^^  and  which  are  0(log(r))  as  r-^.  However,  since  p^^^  must 

decay  at  least  as  fast  as  which  is  0(l/r)  as  r->0D,  we  see  that  this 

cannot  be  allowed.  Consequently,  we  shall  now  require  that  D^p^^^  and  D^u^^^ 

2 

are  0(l/r  )  as  r-w,  and  we  will  show  that  this  assumption  allows  us  to 
determine  the  quantities  of  interest  in  a  consistent  manner.  Using  (7.2),  we 
can  express  this  assumption  as 


D^p^^^  =  0(l/r^) 


=  (l/r^)[rD^f^"'^^^]l^(cos(5))cos(m\('-/?)  +0(l/r^), 


where  is  bounded  as  r-^.  Also,  as  we  shall  show  below,  it  is 

convenient  to  think  of  the  argument  t’  of  f  as  being  replaced  by  a  new 
argument,  say,  r  =  r(t-r+l,tf,^,rj^,r2, . . .  ,M) ,  which  must  be  determined.  Then 

we  can  write 

(7.4a)  Djf^"'^^^(r)  =  f (D^r) . 
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Furthermore,  the  angular  dependence  of  the  terms  on  the  right  sides  of 


equations  (7.3)  suggests  that  D^r  should  be  proportional  to  P™{cos(0))' 
cos(m^-/?).  Thus,  we  write 


(7.4b) 


=  h(r)P“(cos(0))cos(m^-/3) , 


(7.5) 


where  h(r)  must  be  determined.  With  these  assumptions,  we  see  that  all  of  the 
terms  on  the  right  side  of  (7.3)  are  0(l/r  )  (at  least)  as  r-»a,  as  required. 
Thus,  we  can  express  the  (particular)  solutions  for 
in  the  form 


(1/r) log(r)FQ  -  (l/r^)log(r)FQ  +  (l/r^)Fj^ 


2 

■  ^F^(  cos(^)cos(m\{'-/3)j  +0(l/r^), 

•  (l/r)log(r)F^  +  (l/r^lFgj 

2 

■  ^P^(co5(0)cos(m^-/?)  j  +0(l/r^), 

•  (l/r^)log(r)FQ  +  (l/r^lFgj 


■2I^(cos(0))p”’ (cos(^))sin(0)cos^(m^-/3)  +  Od/r^*), 


(l/r^sin(^))iog(r)FQ  +  (l/r^sin(fl) )F^ 


2 

■2m|^P^(cos(^) )j  cos(mf-/?)sin(m^-/9)  +0(l/r'^). 


where  each  F.  =  F.(r,r)  is  a  function  to  be  determined  and  the  primes  denote 
3  3 

differentiation  with  respect  to  r.  Substituting  the  expressions  (7.5)  into 
equations  (7.3),  and  using  (7.4),  we  obtain  the  relations: 


(7.6a)  -  rf-^^^^h(r;  +  0(l/rj, 

(7.6b)  =  ^2  ((3-7)/2)f^”'"^^f^"^^^  +  0(l/r), 

2 

(7.6c)  =  Fg  +  ((7-3'y)/2)[f^"'^^^|  +  0(l/r), 

2 

(7.6d)  F^  =  F^  +  ((7-'y)/2)|^f^"'"^^]  +0(l/r). 

In  (7.6),  the  function  F2  can  be  determined  by  examining  the  terms  in 

3 

equations  (7.3)  that  are  0(l/r  ).  From  the  form  of  the  solution  (7.5),  we  see 

(2) 

that  we  must  demand  that  F^  =  0(l/log(r)),  as  r->flD,  since  otherwise  p  will 

not  decay  as  fast  as  the  Od/r)  rate  of  decay  of  p^^\  Using  this  condition, 
we  see  that  the  right  side  of  equation  (7.6a)  must  be  0{l/log(r)),  as  r->®. 
Using  the  relation  (7.4b),  this  condition  can  be  written  as 

(7.7)  D^(g)  +  (('K+l)/4)P^(cos(^))cos(m#-/3)(l/r^)DQ(r^g^)  =  0(  1  / ( r^log( r ) ) ) , 

as  r^,  where  g  =  f^^^^^/r.  (Here  we  have  assumed  that  D^r  =  -1  + 

0(l/log(r))  and  drldt  =  1  +  0(l/log(r)),  as  r->oo.  See  equations  (7.9)  below.) 
Equation  (7.7)  is  the  major  result  of  this  section  and  it  serves  to 

determine  the  behavior  of  functions  of  r ^ .  It 

can  be  interpreted  as  a  "spherical"  Burgers’  equation  in  the  "time-like" 
variable  ((7+l)/2)P^(cos((?))cos(m^-^)rj.  In  particular,  we  can  express  the 

solution  of  (7.7)  which  satisfies  the  "initial”  condition  g  =  f ^ { t+l-r ) /r 
when  M=0  as 

(7.8a)  g  =  f^"'^^^(T)/r,  where 
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(7.8b) 


r  =  t-r+1  +  ( ('y+l ) /2)F^(t ,  r  j  )P^(  cos(  0) ) cos ( , 
n+1 


=  -V' ^  ^l/s)jds 

n  1  (L— 1  n,j  Jy 


=  Mf 


J=1 
(n+1 ) 


n+1 


(T)log(r  /M)  +  '  a  ^  ^  ( r )  I  (M/ r  ) ■^ 

*  4mw4  ^  »  J  N 

j=2 


where  j  =  2,  3,  ...  ,  n+1. 


Thus,  both  the  and  behavior  of  completely 

determined  by  equations  (6.9),  (6.1),  with  t’  replaced  by  t,  and  (7.8b). 

Before  we  consider  some  specific  applications  of  these  results,  we  note 
that,  from  its  definition  in  (7.8b),  r  has  the  following  properties; 


(7.9a)  T  t-r+1,  as  M-K); 

(7.9b)  r  -»  t-r+1,  as  r-*l; 

(7.9c)  drldt  =  1  +  O(l/log(r)),  drldr  =  -1  +  0(l/log(r)),  as  r-+0D; 

{7.9d)  dr/dS  =  0(M),  and  dT/d4'  =  0(M),  as  M-K); 
and 

(7.9e)  Dj^r  =  ((7+l)/2)f^(cos(tf))cos(ni^-/9lf^'^^^^(r)/r  +  0(l/r^),  as  r-+®. 


Properties  (7.9a)  and  (7.9b)  insure  that  the  linear  solution  is  recovered, 
either  as  SHO  or  as  we  approach  the  surface  of  the  sphere.  Property  (7.9c) 
shows  that  our  definition  of  r  is  consistent  with  the  assumptions  made  above, 
while  property  (7.9d)  shows  that  the  angular  variation  of  r  is  ''small"  and, 
hence,  will  be  described  by  the  next  level  of  the  perturbation  expansion. 
Property  (7.9e)  shows  that  is  0(l/r),  as  r-^,  and  hence  r  satisfies  the 

requirement  that  rD^r  is  bounded  as  r^.  We  shall  comment  further  on  the 

properties  of  t  in  section  9. 
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8. 


As  our  first  application,  we  consider  the  pulsating  sphere  (see  Figure 
2(b)).  In  this  case,  the  (nondimensional )  radial  component  of  the  velocity  of 
the  surface  of  the  sphere  is  given  by  V  =  sin(wt)  and  hence  we  set  n=0  and  m=0 
in  our  results  from  sections  6  and  7.  Then  the  steady  state  solution  for  f(t) 
is  given  by  (6.10)  with  n=0,  and  (p  is  given  by  (6.9)  with  n=0,  i.e. 

(8.1)  f(t)  =  (1+0)^)  ^l^sin(iJt)  -  c<;cos(wt)j  and  p  =  -f(r)/r, 
where  r  is  determined  by  (7.8b)  with  n=m=0,  i.e., 

(8.2)  T  =  t-r+1  +  ((7+l)/2)Mf ’(r)log(r^/M) . 

Using  these  definitions  in  equations  (6.1)  and  (7.1)  we  find  for  this  case 
=  (^w/((l+w^)r)j  (cos(u;r)  +wsin(«r)), 

(2) 

p  =  2bgWcos(2wr)/r  -  2b2wsin(2u;r)/r 

-  l/(4r‘^(l+<*;^) )  +  (l-7)u;^/(4r^(l+u»^)) 

+  ((17+7)w^/(l6r^(l+w^)^)  -  ( 13-37) (w^-l)/(64r'*(l+a;^)^) 

(8.3)  +  (5+7)u;^(w^-l )  /  (8r^(  l+w^)^) j  cos(2wt ) 

+  [(13-37)w/(32r'*{l+w^)^)  -  (5+-)w^/(4r^f  1+w^)^) 

+  (17+7)a;(w^-l )  /  (32r^(  l+w^)^)j  s  in(2u;r) 

+  0(l/r®), 
and 

bg  =  -3a;^(14-27+17w^+7w^)/{16(l+w^)^(l+4w^)), 
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bg  =  w(29-3"?-17u;*^  +  15iw^-64u;'*)  /  (32S  1 


»  j 


(The  corresponding  expressions  for  the  perturbation  coefficient  functions 
and  for  the  radial  component  of  the  velocity  are  listed  in  the 

appendix .  ) 

In  Figure  3,  we  have  plotted  r(p^^*+  isolid  lines)  us  a  function 

of  r  at  a  fixed  time  for  iv=1.25  and  for  two  different  values  of  M  Here  it  is 

convenient  to  think  of  '  as  our  approximation  to  the  norma i i zed , 

nondimens ional  acoustical  density  field  (p-l)/M  (from  equation  (3  2)i.  In 
Figure  3  we  have  also  plotted  the  approximation  to  r(p-l)/M  based  on  Whitham’s 
tl974)  first  order  solution  (long  dashed  lines),  the  classical  linear  solution 
(short  dashed  lines),  and  a  solution  obtained  by  purely  numerical  methods 
(circles).  This  numerical  solution  was  obtained  using  a  MacCormack 
predictor-corrector  scheme  oii  equations  (5.1)-{5,4).  As  Figure  3  illustrates, 
there  is  very  good  agreement  between  the  multiple  scales  perturbation 
solutions  and  the  numerical  solutions,  even  for  a  Mach  number  of  0.3,  where 
the  gradual  steepening  and  asymmetry  of  the  wave  is  apparent, 

b)  The  Oscillating  Sphere 

As  a  second  application,  we  consider  the  oscillating  sphere  (see  Figure 
2(c).)  In  this  case,  the  radial  component  of  the  velocity  of  the  surface  of 
the  sphere  is  given  by  V  =  cos(^)sin(wt)  and,  hence,  we  set  n=l  and  m=0  in  our 
results  from  sections  6  and  7.  Then  the  steady  state  solution  for  f(t)  is 
given  by  (6.10)  with  n=l ,  and  p  is  given  by  (6.9)  with  n=l ,  i.e. 

(8.4)  f(t)  =  (4+1*;“*)  ^  |^(2-w^)sin(wt)  -  2wcos(u;t)j  , 

p  =  -  cos(fl)  l^f  ’  (r)/r  +  f(7-)/r^j, 

where  r  is  now  determined  by  (7.8b)  with  n=l  and  m=0,  i.e. 
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(8.5)  T  =  t-r+1  +  ((7+l)/2)M|f”(r)log(r^/M)  -  f  ’  ( r )  j^(M/ r  ^ ) -1 J  |cos  ( 0)  . 

V  / 

Using  these  expressions  in  equations  (6.1)  and  (7,1),  and  using  Wolfram’s 
(1991)  symbolic  computation  system  Uathemat  ica ,  we  find  for  this  case 

^u(2-u;^+2ru^)cos(uir) /r^  +  2-2r+ru;^) s inlt^r ) /r^j  eosCfl )  /  (4+^"* ) , 

(2)  f  q 

p  =  j2a;(fQ^cos(2cjr)-f^^sin(2wr))/r  -  {3w/2)(  f2gCos(2wT)+f2^sin(2u>r)  )/r 

+  f2pCos( 2£jr  )+f2gSin(2cjr) ) /r""  +  2w^(  f^^cos (2wt l-f^^s  i n(2cjr) ) /r 

+  /Sr^  +  ( l-'7)w^(4+w'^)/8r^ 

+  [(323+37)w^(2-<j^)/64r®  -  (7+'r)w'‘{2-w^)/2r^ 

-  (259+19'7)a;^(4-8a;^-Hi;*^)/192r^  +  (5+'7)«'^(4-8w^+w"^)/16r^jcos(2u;T) 

+  (^(5+'y)w®(2-(j^)/4r^  -  (259+19'y)u;^{2-w^) /48r'‘ 

+  (7+7)w^(4-8w^+w‘*)/8r^ 

-  w(  323+37)  ( 4-8w^+u;'^ )  / 256r®]  s  i n(  2wt  ) | /  ( 4+w ' ) ^ 

+  cos(20)|(9w/2)  (-f2gCos(2wr)+f2^sin(2wr') )  /  r^ 

-  9w^(f2^cos(2wr)+f2gSin(2wr))/r^  +  Sw^C f2gCOs(2wr)-f2^sin(2wr) ) /r 

+  (3-7)u;^(4+w^) /8r^  +  ( (l-7)w‘*(4+w^) /8r^ 

+  [3(33+7)w^(2-w^)/32r®  -  (17+7)w^(4-8a;^+u;'^)yi6r'^ 

-  3(9+7)w^(2-a;^) /8r^  +  (5+7)w^(4-8w^+4i!^) /16r^j  cos(2wr) 

+  ^(5+7)w^(2-w^) /4r^  - 

+  3(9+7)w^(4-8w^+w^)/32r^  + 
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(2) 

The  constants  f2g,  and  which  appear  in  the  definition  of  p 

are  related  to  the  homogeneous  solution  to  equations  (7.1)  and  are  defined  in 

the  Appendix.  (The  form  of  the  solutions  for  u^^\  ul^\  and  ul^^are 

r  p  r  (f 

also  presented  in  the  Appendix.) 

The  form  of  these  expressions  allows  us  to  make  an  interesting 
observation  concerning  the  angular  dependence  of  the  sound  radiated  from  an 
oscillating  sphere.  Specifically,  it  is  interesting  to  examine  the  density 

and  velocity  in  the  d=7r/2  direction.  In  this  direction,  both  a* d 

are  zero,  while  is  nonzero.  Jn  contrast,  for  this  direction  both  p^^^ 

and  u^^^  are  nonzero,  while  ul^^  is  zero  (since  it  contains  a  factor  of 
r  u 


2 

sin(2d)).  Thus,  there  is  a  ’'small"  (0(M  ))  amount  of  sound  radiated  in  this 
direction,  which  has  a  characteristic  frequency  of  2w.  In  particular,  in  the 
far  field  we  can  use  the  expressions  above  to  write 


p^2^=  uj.^^  Od/r^) 


=  (2w/r)  |^(2f2gW^-fQg)cos(2u;r)  +  ( fQp-2f2j.w^) sin(2uir) j  +  0(l/r^), 


along  $=itl2. 


In  addition,  we  find  that  the  expressions  for  and  correspond  to  a 

rotational  flow,  even  though  u|^^  and  u^^^  correspond  to  an  irrotat  ional  flow. 


The  presence  of  vorticity  in  the  flow  can  be  interpreted  as  being  due  to 
tangential  acceleration  of  the  fluid  at  the  boundary,  as  discussed  by  Morton 
(1984).  This  vorticity  appears  only  very  near  the  sphere  and  is  given  by 
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Vxu  =  i^Msin(2d)||^(7+l)w'*/(4(4+t«;'^))j  ( log(r)  )/r^ 

+  4-8w^+ti;’^ )  /  ( 4+w'^ ) cos ( 2wr )  / 

+  ^4«^(2-a;^) /(4+w^)^j  s  in(2wr )  /  +  Odyr"*) 


In  Figure  4,  we  have  plotted  approximations  to  r{p-l ) / (Mcos( 0) )  as  a 

2 

function  of  r  at  a  fixed  time  for  0=0  and  0=vl3,  as  well  as  np-ll/M  at 
0=ir/2,  with  w=1.25  and  for  two  different  values  of  M.  These  approximations 

are  the  multiple  scales  solution  r(/j^^^+Mp^^^ )/cos(d)  (solid  lines),  the 
classical  linear  solution  (short  dashed  lines),  and  a  solution  obtained  by 
purely  numerical  means  (circles).  For  this  case,  we  found  that  numerical 
approximations  based  on  a  Jameson  type  Runge-Kutta  finite  volume  scheme, 
which  is  second  order  accurate  in  both  space  and  time,  gave  better  agreement 
with  the  multiple  scales  solution  than  numerical  approximations  based  on  the 
MacCormack  scheme,  as  in  the  case  of  the  pulsating  sphere.  (We  thank  Mr. 

David  Lockard  and  Dr.  Kenneth  Brentner  of  the  NASA  Langley  Research  Center  for 
carrying  out  these  calculations  for  us.)  Again,  the  agreement  between  the 
multiple  scales  solutions  and  the  numerical  solutions  is  qualitatively  good. 

c)  A  "Squishing"  Sphere 

As  a  third  example,  we  consider  the  case  when  the  radial  component  of  the 

2 

surface  velocity  of  the  sphere  is  given  by  V  =  sin  (0)sin(2^)sin(u;t)  (see 
Figure  2(d)).  Thus,  we  set  n=2  and  m=2  in  our  results  of  sections  6  and  7, 
and  find  from  (6.10)  with  n=2  that  the  steady  solution  for  f(t)  is  given  by 

(8.7)  f(t)  =  (3(81+9w^-2w^+w^) )  ^  ^(9-4w^)sin(wt)  -  w(9-<*(^)cos(wt)j  , 
while  the  solution  for  ip  is  given  by  (6.9)  with  n=2,  i.e., 

(8.8)  ip  =  -[f”(r)/r  +  3f’{r)/r^  +  3f(r)/r^]3sin^(5)sin(2^) , 
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with  T  determined  by  (7.8b)  with  n=2,  i.e.. 


T  -  t-r+l  +  (('7+l)/2)F2(r,r^)3sin  (tf))sin(2^)  , 
.(T.r^)  =  M|f” '(r)log(r^/M)  -  3f  ”  (r)  f(M/r^)-ll 

-  (3/2)f’(r)[(M/r^)^-l 


Using  these  expressions  we  find 
=  sin(2^)sin^(5)  1^ 

( (27tj-12a;^)/r^+(27w^-3w^)/r^-(9(j^-4w^)/r)cos(u>T)  +■ 

( (27u^-3w'^)  /r^- (27cj^-12u;*^) /r^-Ow'^-w^)  /  r ) sin(«T )  j  /(81+9a;^-2u^-KJ®) , 

while  expressions  for 

order  perturbation  coefficient  functions,  are  given  in  the  Appendix. 

2 

In  Figure  5  we  have  plotted  approximations  to  r(p-l)/(Msin  (0)sin(2^)) 
as  a  function  of  r  at  a  fixed  value  of  time  with  M=0.6  for  =  x/12,  w/6,  and 
7r/4  at  (a)  5=w/4  and  (b)  8='k/2,  using  the  multiple  scales  solution 
(1)  -  2 

Tp  /(sin  (^)sin(2^))  (solid  lines)  and  the  classical  linear  solution 
(dashed  lines).  Due  to  the  excessive  storage  requirements  for  this  fully 
three-dimensional,  time  dependent  problem,  no  numerical  approximations  were 
computed. 

9.  Conclusions  and  Discussion 

In  this  section  we  shall  first  make  some  general  observations  concerning 
the  perturbation  approach  we  have  employed  to  the  problem  of  sound  generation 
by  vibrating  bodies,  and  then  comment  specifically  on  the  insights  gained  from 
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the  examples  we  have  presented. 

We  note  first  that  the  multiple  scales  perturbation  approach  has  allowed 
us  to  "capture"  many  of  the  salient  nonlinear  features  of  acoustic  wave 
motion.  For  example,  as  the  Mach  number  increases,  the  nonlinear 
characteristics  of  the  gradual  steepening  of  the  waves,  the  asymmetries  due  to 
"second  harmonic"  terms,  as  well  as  the  convergence  of  the  wave  to  an  "N-wave" 
profile,  are  all  evident  in  the  examples  we  have  considered.  These  phenomena, 
of  course,  are  not  predicted  by  linear  theory.  Thus,  it  appears  that  the 
method  has  allowed  us  to  obtain  approximate  analytical  solutions  which  are 
valid  for  (much)  larger  values  of  M  than  are  the  classical  linear  acoustic 
solutions.  In  addition,  the  solutions  assume  a  rather  simple  form,  being 
essentially  of  the  same  form  as  the  classical  linear  solutions,  but  with  a 
different  argument.  That  is,  the  classical  retarded  time  t+l-r  in  the  linear 
solution  has  simply  been  replaced  by  r,  where  r  is  defined  implicitly.  Other 
perturbation  approaches  to  this  problem  have  resulted  in  more  involved 
expressions  for  the  quantities  of  physical  interest  (cp.,  e.g..  Crow  (1970)). 
(We  should  note  that,  if  the  present  method  of  analysis  were  to  be  continued 
to  a  higher  order  than  that  presented  here,  various  terms  of  the  form 

M^|^log(M)j  should  be  added  to  the  expansions  (5.7).  However,  since  the 

2 

present  analysis  was  terminated  at  terms  which  are  0(M  ),  these  terms  were  not 
needed  and,  hence,  were  not  shown  explicitly.) 

The  general  form  of  the  equation  for  t  also  provides  an  insight  into  the 
interplay  between  the  effect  of  the  Mach  number  and  the  angular  variation  of 
the  acoustical  source  on  the  form  of  the  acoustic  wave.  For  example,  if  M=0 

or  if  the  term  P^(cos(0) )cos(m^-/9)  is  zero  (corresponding,  say,  to  a 

particular  value  of  ^  or  ^),  then  from  (7.8b)  we  see  that  T=t-r+l  and,  hence, 
our  solution  (to  leading  order)  reduces  just  to  the  linear  acoustic  solution. 

However,  if  M?^0,  but  the  angular  terra  P^(cos(0)  )cos(m^-/?)  is  zero,  then  our 

solution  (to  leading  order)  again  reduces  to  the  linear  solution  and  hence  the 

nonzero  Mach  number  has  only  a  second  order  (i.e.  0(M  ))  effect  on  the  form  of 
the  solution.  These  observations  motivate  us  to  define  an  "effective"  Mach 
number  M^^^  by 
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M  =  MP™(cos(0))cos(in\i'-/3)  . 
et  I  .  n 

Thus,  it  is  the  (angularly  dependent)  effective  Mach  number  ,  and  not 

just  M,  which  determines  the  magnitude  of  departure  of  the  wave  from  its 
linear  form.  This  phenomena  is  illustrated  in  Figure  4  for  the  example  of  the 
oscillating  sphere,  for  which  n=l  and  m=0  in  the  formulas  above,  and  also  in 
Figure  5,  corresponding  to  n=2  and  m=2.  In  Figure  4(a),  for  example,  the 
linear  form  of  the  wave  is  clearly  seen,  even  for  nonzero  Mach  numbers,  when 
cos(5)  approaches  a  value  such  that  is  zero. 

In  comparing  our  results  with  Whitham’s  (1974)  results,  we  note  that  for 
spherically  S3anmetric  disturbances  (for  which  n=0)  our  equations  (7.8)  reduce 
to  his  results.  However,  in  more  general  problems  (for  which  n>0),  our 
results  differ  from  Whitham’s  equations.  In  the  far  field,  where  the  log  term 

in  the  expression  for  F^  dominates  the  remaining  terms,  equation  (7.8b) 

reduces  to 

r  =  t-r+1  +  ((7+l)/2)f^'’'''^^(r)log(rj^/M)MP^(cos(0))cos(m^-/?) . 

For  n=0,  this  equation  is  equivalent  to  Whitham’s  equation  for  t.  For  n>0, 
we  use  r^^/M  =  r  to  note  that  this  equation  is  the  same  as  Whitham’s  result 

when  we  replace  M  in  his  formula  by  . 

For  n>0,  the  definition  of  r  involves  terms,  through  the  definition  of 

F^,  which  are  negligibly  small  in  the  far  field.  However,  these  terms 

appear  to  have  an  interesting  (and  vital)  interpretation.  In  particular,  the 

terms  other  than  the  log  term  in  F  give  rise  to  terms  in  D, r  which  are 

"  n  "  1 

2 

0(l/r  )  (at  least)  as  r-»oo,  which,  in  turn,  give  rise  to  terms  on  the  right 

3 

sides  of  equations  (7.1a,b)  which  are  0(l/r  ),  as  r-^.  However,  without  these 
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terms  present,  no  steady  state,  periodic  (pariicular)  solution  exists  for 

p^^\  this,  we  note  first  that,  if  a  periodic 

solution  with  period,  say,  T,  for  these  quantities  does  exist,  we  can 
integrate  equations  (7.1b)-{7. Id)  over  a  complete  period  in  t  and,  with  a 
little  manipulation,  obtain  the  relations 


0  =  idldB) 


[fright  si 


side  of  (7.1b)  dt  -  D 


0, 


fright  side  of  (7.1c)jdt, 


oT, 


0  =  0/5^)J  right  side  of  (7.1b)jdt  -  j^right  side  of  (7.1d)|dt 


The  first  term  on  the  right  side  of  each  of  these  relations  involves 

while  the  second  terms  in  each  relation  is  independent  of  For 

example,  in  the  case  of  the  oscillating  sphere,  the  first  of  these  relations 
yields  the  requirement 


0  =  (aj^  2  l)('lf+l)w^/(8r^(4+y'*) ) . 


Thus,  if  a  periodic  solution  is  to  exist,  we  must  have  a  „  =  -1  =  ^,  as 

X  f  ^  1  f  ^ 

stated  in  equation  (7.8b).  Similar  equations  hold  for  n^2,  which  led  to  the 
particular  definition  of  which  appears  in  equation  (7.8b). 

It  is  also  interesting  to  note  how  the  distortion  of  the  acoustical  wave 
from  a  linear  wave  varies  with  the  parameter  n.  In  Figure  6  we  have  plotted 

r(p-l)/M  at  M=0.3  for  (a)  the  pulsating  sphere  (n=0),  (b)  the 
oscillating  sphere  at  ^=0  (n=l),  and  (c)  the  squishing  sphere  at  $-Ttl2  and 
^=y(IA  (n=2).  (In  (b)  and  (c)  the  angular  variables  were  chosen  so  that  the 
distortion  of  the  wave  was  maximum.)  As  the  figure  illustrates,  for  a  given 
value  of  M  the  amount  of  distortion  decreases  as  n  increases.  Intuitively, 
we  may  think  of  this  as  being  due  to  the  increasing  number  of  "degrees  of 
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freedom”  of  the  wave.  Alternatively,  to  obtain  the  same  amount  of  distortion 
in  the  wave  as  n  increases,  the  value  of  M  would  also  have  to  be  increased. 
This  observation  is  consistent  with  our  one-dimensional  example  (which  can  be 
thought  of  as  a  very  "confining"  flow),  since  for  this  example  there  was 
significant  distortion  of  the  wave  at  a  much  lower  value  of  the  Mach  number. 

We  also  note  that  the  analysis  presented  here  may  provide  some  insight 
into  the  derivation  of  appropi iate  non-reflecting  boundary  conditions  to  be 
used  at  artificial  computational  boundaries  for  a  purely  numerical  simulation 
of  acoustical  waves  generated  by  general  (arbitrary)  sources.  For  example, 
the  rather  simple  form  of  our  final  results  allows  us  to  express  the  far  field 
behavior  of  the  waves  in  a  concise  form.  This  expression,  when  used  with  some 
of  the  ideas  of  Bayliss  and  Turkel  (1980),  for  example,  may  allow  us  to  derive 
the  required  boundary  conditions.  Investigations  along  these  lines  are 
continuing. 
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Appendix 

In  this  appendix  we  record  the  form  of  the  perturbation  coefficient 
functions  for  the  velocity  components,  as  well  as  various  constants, 
associated  with  the  examples  presented  in  section  8. 

Pulsating  sphere: 

ujO)  ^  j^(r-l)w/((l+a;^)r^)]cos(wT)  +  (l/(l+u;^))  .;^/r  +  1  / r^j  s in(wr) , 

=  (b_/r^  +  2b_w/r)cos(2wr)  +  (b  /r^  -  2b„w/r)sin(2ur) 
r  A  o  o  ^ 

+  [(29-3'y)w^/(l6r^(l+w^)^)  +  w^(a;^-l)/(r^(l+w^)^)]  cos(2w7-) 

-  ((l+7)a;^log(r)  )/(4r^(l-K(;^)) 

+  ^(29-3'y)a;(a;^-l}/(32r^(l+u;^)^)  -  2u;^/(r^(l-HJ^)^)]  sin(2wr) 

+  0(l/r^) . 

Oscillating  sphere: 

|^2(r-l)w(2+ra;^)cos(wr)/r^  +  (4-2w^+4rw^-2r^w^+r^u;^)sin(wr) /r^j 

•  cos(0)/(4+w'^) , 

^(w(2r-2-rw^)cos(ci;T) /r^  +  (2-w^+2rci;^) sin(wr) /r^j  sin(5)  / (4+^"^) , 

KQ(r)+R^(r)cos(2u^)+R2(r)sin(2c.;r)  + 

^SQ(r)+Sj^(r)cos(2u;r)+S2(r)sin(2w7)j  cos (2^) , 

^TQ(r)+Tj(r)cos(2wr)+T2(r)sin(2wT)j  sin(2fl) , 
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fQ^=  [^2120+727+(  1106-430'y)a;^+(519+7'7)w^+(408+24'7)cj^  j  /  ( 384(  1+4[J^)  )  , 

fog=  -a;[’2028-207+(344-6807)w^-(l037-4997)w'^ 

+(360+1047)w®-512u;^j  /(768(l+4w^)  )  , 

f2^=  w^|^-2106+1987-(141-5637)a)^+(407-5537)w'‘-(280+1047)w® 

+  (408+247)!J®j  /(24(81+36tJ^-32u;‘^+64w®) )  , 

f2g=  wj^4140-4687-(6136+2487)u;^-(4261-329l7)w'^+ 

(96-4807)a;®-(2512+3367)a;®+512w^^j  /I96(81+36u;^-32u;^+64w®) ) . 

Here  each  R. ,  S.  and  T.  can  be  expressed  as  a  multiple  of  (log(r))/r  and  a 
J  J  J 

power  series  in  the  variable  (1/r),  beginning  with  a  term  proportional  to  1/r 

Squishing  sphere: 

u^^^=  sin(2^)sin^(^)  (4w^-9u;^) /r  +(36w^-4(*;^)/r^ 

+  (81w-36u;^)  /r^-(81w-9w^)  /r'^lcoslwt)  + 

( { 81-36w^ )  /  r ( 8lu;^-9w'* )  /  r ^-  ( )  /  r ^-  ( 9w‘*-u;® )  /  r )  s  i  n ( wr ) )  j 

/(81+9u;^-2c.;^-kj^)  , 

u^^^=  2sin(2f)sin(^)cos(d) 

( (27u)-3<*;^)/r'^-(27w-12w^)  /r^-(9u;^-w^)/r^)cos(wr)  + 

( ( -27+12w^ )  /  r  ( 9w^-4w^ )  /  r^-  ( 27w^-3w^  )/r^)sin(wr)j/(  81+9a;^-2w'^+a;® ) , 
u^^^=  2cos(2^)sin(5) ^ 
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4- 


((27u;-3«^)/r'^-(27u-12u^)/r^-(9w^-tj^)/r^)cos(wr) 
((-27+12w^)/r^+(9w^-4w'‘)/r^-(27w^-3w'^)/r^)sin{t^r)]  /(8U9w^-2u;'‘+w^j . 


2  .  1 

^  'cos(2k^)  H  cos(4ji/')  ■‘'’^k  j‘^°®^2wr)+RSj^  .sin(2wr) 
k=0  '^j=0 


2/1  . 
u|.^^=  ^  cos(2k^)p  '  cos{4j^)^Sj^  j'^'^k  j*^®®^2wr)+SSj^  .  s  in{ 2u;t' )]  | , 
k=0  ^j=0  ’  ’  ’  '' 


2/1  s 

Ug^^=  ^  ^sin(2k^)n  cos(4j^)[T^^  .cos(2t4;r)+TSj^  .  s  in(2u;r)j  i , 

k=l  ''j=0  ’  ’  ’  ' 


^^sin((2k-l)0)|sin(4t/')  l'''^^k  i‘^°®^2wr)+USy^  j^sin(2ur)j 

k=l  I  ’  ’ 

+  Uj^  0^^^  qCos(2wt)+US^  Qsin(2wT) 


Here  each  of  the  coefficients  R,  RC,  .,  ...  ,  US,  .  can  be  expressed  as  a 

J  j  K,  j 

multiple  of  (log(r))/r  and  a  power  series  in  the  variable  (1/r),  beginning 
with  a  term  proportional  to  1/r. 
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fjjgur£_^:  Approximations  to  (/?-l)/M  for  the  simple  one-dimensional 

example  plotted  as  a  function  of  distance  from  the  source  (at  x=0)  for 
t=42.44,  using  the  multiple  scales  approximation  from  equation  (4.18)  (solid 
lines),  the  classical  linear  solution  (4.15)  (short  dashed  lines),  and  a 
solution  to  equations  (4.1)-(4.2)  obtained  by  purely  numerical  means 
(circles),  with  (a)  M=0.01  and  (b)  M=0,04. 


Figure  2:  (a)  An  illustration  of  the  Cartesian  and  spherical  coordinate 

systems  used  in  sections  5-8.  Also,  an  indication  of  the  vibrations  of  the 
sphere  for  the  three  examples  considered  in  section  8:  (b)  the  pulsating 
sphere;  (c)  the  oscillating  sphere;  and  (d)  the  "squishing"  sphere,  indicating 
the  motion  of  the  sphere  in  the  plane  ^=7r/2. 
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Squishing  S 


w=l .25 

t=50 

Squishing  Sphere" 

M=( 

Figure  5:  Approximations  to  r(p-l)/(Msin  {6)&in{2f))  for  the  "squishing 

sphere"  example  plotted  as  a  function  of  distance  from  the  center  of  the 
sphere  at  t=50  with  u=l.25  and  M=0.6,  for  (a)  0=k/4  and  (b)  B=n/2,  at  ^  = 

?r/12,  7r/6,  and  ir/4,  using  the  multiple  scales  solution  rp  /(3sin  (0)sin(2^) 
(solid  lines)  and  the  classical  linear  solution  (dashed  lines). 
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